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Practical aspects of identifying dynamic models for aircraft in real time were studied. 
Topics include formulation of an equation-error method in the frequency domain to estimate 
non-dimensional stability and control derivatives in real time, data information content for 
accurate modeling results, and data information management techniques such as data 
forgetting, incorporating prior information, and optimized excitation. Real-time dynamic 
modeling was applied to simulation data and flight test data from a modified F-15B fighter 
aircraft, and to operational flight data from a subscale jet transport aircraft. Estimated 
parameter standard errors and comparisons with results from a batch output-error method 
in the time domain were used to demonstrate the accuracy of the identified real-time models. 
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body-axis translational accelerometer measurements, g or ft/sec - 
wing span, ft 

mean aerodynamic chord, ft 

body-axis non-dimensional aerodynamic force coefficients 
body-axis non-dimensional aerodynamic moment coefficients 
expectation operator 
mass moments of inertia, slug-ft 2 

imaginary number = x[-\ 
cost function 
aircraft mass, slug 

body-axis pitching moment from engine thrust, ft-lbf 

body-axis roll, pitch, and yaw rates, rad/sec or deg/sec 

dynamic pressure, lbf/ft 2 

standard error 

wing reference area, ft 2 

maneuver length, sec 

body-axis engine thrust, lbf 

airspeed, ft/sec 

angle of attack, rad or deg 

sideslip angle, rad or deg 

elevator, aileron, rudder, and trailing-edge flap deflections, rad or deg 

canard, differential canard, stabilator, and differential stabilator deflections, rad or deg 

Euler roll, pitch, and yaw angles, rad or deg 
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0 

= parameter vector 

27 

= covariance matrix 

CO 

= frequency, rad/sec 

superscripts 

T 

= transpose 

f 

= complex conjugate transpose 

A 

= estimate 


= time derivative 

-1 

= matrix inverse 

subscripts 

o 

= reference value 


I. Introduction 

D YNAMIC modeling in real time has many important practical uses, such as improving the efficiency of 
stability and control flight testing, flight envelope expansion, adaptive or reconfigurable control, vehicle health 
monitoring, and fault detection. Several methods 1 5 have shown potential for identifying local linear dynamic 
models from flight data in real time. One of these methods 46 is based on a recursive Fourier transform and 
equation-error modeling in the frequency domain. This method, sometimes called the Fourier Transform Regression 
(FTR) method, produces very accurate results with valid error measures and has many practical advantages. The 
FTR method has also been independently evaluated 7 ' 8 as the best method available for real-time dynamic modeling. 
For these reasons, the FTR method was chosen for further study and application. 

The FTR method has been successfully applied 611 to identify accurate linear dynamic models in real time at 
individual flight conditions. While this capability is important and useful, further progress requires that this 
capability be extended to continuous application as the aircraft flies through a wide range of changing flight 
conditions throughout the flight envelope. Ultimately, local real-time modeling results could be integrated into a 
global aerodynamic model that could be updated in real time as the aircraft changes flight conditions, changes 
configuration, ages, or becomes damaged in some way. This vision of real-time global dynamic modeling has many 
important implications for efficient flight testing, accurate flight simulation, adaptive or reconfigurable control, and 
aircraft safety. 

One important aspect of applying real-time dynamic modeling for varying flight conditions and aircraft 
configurations is determining the data information content requirements for accurate dynamic modeling results. 
Changing aircraft flight conditions or aircraft configurations means that parameters in the approximating dynamic 
model change. Dynamic motion of the aircraft, either from ordinary flight operations or from control surface 
excitation, is necessary so that the measured data will exhibit the aircraft dynamics to be modeled. Naturally, if the 
real-time dynamic modeling is to be done continuously or on a regular basis, it is important that only the minimum 
necessary aircraft excitation be applied, and the resulting aircraft motion should be as small and unobtrusive as 
possible. 

This paper investigates data information requirements for accurate real-time dynamic modeling. Flight 
experiments on a modified F-15B fighter aircraft are used to illuminate issues related to data information content 
necessary for accurate real-time modeling. Real-time modeling is also applied to operational flight data from a 
subscale jet transport model, to evaluate the feasibility of real-time modeling without specific excitation. This is an 
important step in extending local real-time modeling to the case of changing conditions, interpreted broadly to 
include flight condition changes, configuration changes, damage, and failure scenarios. Issues such as data 
information content necessary for fast and accurate local modeling, model validation, necessary excitation, data 
forgetting, and methods for incorporating prior information are studied. 

The next section describes the methods used. A model formulation is developed that retains full nonlinear 
dynamics, with linearized aerodynamic models. The FTR method is described, along with explanations of methods 
for data forgetting and incorporating prior information into the real-time parameter estimation algorithm. Next, the 
flight test aircraft are described, including flight instrumentation and characteristics of the flight data. The results 
section includes simulation and flight test investigations examining data information requirements for accurately 
identifying local dynamic models in real time. Finally, the concluding remarks summarize progress made so far and 
outline next steps. 
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II. Methods 


A. Aerodynamic Modeling 

Non-dimensional aerodynamic force and moment coefficients for an aircraft can be computed from flight 
measurements as follows 6 : 
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These expressions retain the full nonlinear dynamics in the aircraft equations of motion. For local real-time 
modeling over a short time period, the force and moment coefficients computed from Eqs. (1) and (2) can be 
modeled using linear expansions in the aircraft states and controls: 
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The A notation indicates perturbation from a reference condition. In each expansion, a single term is shown to 
represent all relevant (and similar) control terms, to simplify the expressions. For example, in Eq. (3c), the term 
C ltl j AS represents terms for each control that affects C m , e.g., C mg AS e + C m AS y . In Eq. (3c), C m ^ 

represents the non-dimensional pitching moment at a reference condition, and similarly for the other expansions. 


3 

American Institute of Aeronautics and Astronautics 



The linear aerodynamic models in Eqs. (3) and (4) contain parameters called stability and control derivatives, 
such as C m and C mg , which characterize the stability and control of the aircraft. For short periods of time, the 

stability and control derivatives are considered to be constant model parameters to be estimated from flight data. 
Repeating the parameter estimation at short time intervals produces piecewise constant estimates for the stability and 
control derivatives, which in general vary with flight condition and changes to the aircraft, such as configuration, 
age, damage, or failures. 

The next subsection describes how the unknown stability and control derivatives in the linear aerodynamic 
models of Eqs. (3) and (4) can be estimated from flight data using equation-error parameter estimation in the 
frequency domain. 

B. Stability and Control Derivative Estimation in the Frequency Domain 

This section describes the FTR method for estimating unknown parameters in a dynamic model in real time. 
Some of the material presented here can also be found in Refs. 4 and 6. 

The first step required for modeling in the frequency domain is to transform the measured flight data from the 
time domain into the frequency domain. The finite Fourier transform is the analytical tool used for this task. For an 
arbitrary scalar signal x ( t J on the time interval [(), 7’] , the finite Fourier transform is defined by 

s -*(®) = | x (t) e * m dt (5) 


which can be approximated by 


jv — 1 
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where x(i\ = x(iAt\ , T = N At , and At is the sampling interval. The summation in Eq. (6) is defined as the 
discrete Fourier transform. 


X(c))=Y J x(i)e- j(aiAt (7) 

i=0 

so that the finite Fourier transform approximation in Eq. (6) can be written as 

x(a>) x X (®) At (8) 

Some fairly straightforward corrections 12 can be made to remove the inaccuracy resulting from the fact that Eq. (8) 
is a simple Euler approximation to the finite Fourier transform of Eq. (5). However, if the sampling rate is much 
higher than the frequencies of interest, as is typically the case for dynamic modeling from flight data, then the 
corrections are small and can be safely ignored. 

The Fourier transform is applied to the non-dimensional force and moment coefficients computed from Eqs. (1) 
and (2) using measured time-domain data. This results in the non-dimensional force and moment coefficients in the 
frequency domain. Often, measurements of the angular accelerations p,q,r are not available. In the frequency 
domain, these derivatives can be calculated easily. For example, the Fourier transform of the rolling moment 
coefficient can be computed as: 
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and similarly for C m and C n . This approach implements the derivative of the body-axis angular momentum in the 
frequency domain, including the time variation in the inertia quantities. Note that the Fourier transform of the 
nonlinear terms is handled by computing the nonlinear terms in the time domain, then applying the Fourier 
transform to the resulting time history. Treatment of the dynamic pressure q in Eq. (9) is consistent with an 
assumption that the dynamic pressure varies slowly, which is a good practical assumption. 

For the perturbation states and controls in Eqs. (3) and (4), time histories of the measured states and controls are 
high-pass filtered to remove the steady part of each signal, then transformed into the frequency domain using the 
discrete Fourier transform. The break frequency for the high-pass filter is set just below the lowest frequency 
selected for the modeling. Fligh-pass filtering is implemented with a fourth-order Butterworth digital filter. 
Similarly, the quantities transformed in Eq. (9) (shown within the square brackets) are also high-pass filtered prior to 
Fourier transformation. This approach effectively drops out the bias terms in the models of Eqs. (3) and (4). The 
high-pass filtering prevents leakage from the relatively large spectral component at zero frequency associated with 
the steady component of the signals from polluting transformed data at low frequencies. 

For each aerodynamic model Eq. (3) and (4), the parameter estimation problem can be formulated as a standard 
least squares regression problem with complex data 6 . 
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where Aq n = IF^qc /(2V and e represents the complex equation error vector in the frequency domain. The 

symbols Aa(k}, k = l, denote the Fourier transform of the angle of attack perturbation state for each 

frequency co k , and similarly for other quantities. Each transformed variable depends on frequency. The 
frequencies co k can be chosen arbitrarily, and are therefore chosen to cover the frequency band where the aircraft 
dynamics lie, as will be discussed later. The least squares cost function is 

J(9)=^{z-xe) f (z-xe) ( 14 ) 
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This cost function contains M squared error terms in summation, corresponding to M frequencies of interest. Similar 
cost expressions can be written for individual lines from Eqs. (3) and (4). The parameter vector estimate that 
minimizes the least squares cost function is computed from 6 
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The estimated parameter covariance matrix is 6 
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where the equation-error variance <j 2 can be estimated from the residuals, 
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and n p is the number of unknown parameters, i.e., the number of elements in parameter vector 9 . Parameter 

standard errors are computed as the square root of the diagonal elements of the C ov^S'j matrix from Eq. (16), using 

<t 2 from Eq. (17). Reference 6 explains why the estimated parameter standard errors are computed in this way, and 
also why this calculation in the frequency domain produces parameter error measures that are consistent with the 
scatter in parameter estimates from repeated maneuvers. Realistic simulation testing has shown that the accuracy of 
model parameters estimated with this method is comparable to using a time-domain output-error method employing 
iterative nonlinear optimization in post-flight batch mode 13 . 

The model formulation given here is widely applicable, because the assumption of constant linear aerodynamic 
models over short time periods is very accurate for non-dimensional stability and control derivatives, where the 
effects of changing dynamic pressure and mass properties are removed. 

To implement this least squares parameter estimation in the frequency domain, the parameter estimation 
calculations in Eqs. (15)-(17) are applied to frequency-domain data at selected times, normally at regular time 
intervals. The frequency-domain data must therefore be available at any time, so the Fourier transforms are 
computed using a recursive Fourier transform, described next. 

C. Recursive Fourier Transform 

For a given frequency a> , the discrete Fourier transform in Eq. (7) at time i At , denoted by Z, (co) , is related to 
the discrete Fourier transform at time (i-\)At by 


X i (co) = X i _ l (co) + x(i)e-j a>iAt 


(18) 


where 


g -ja>iAt _ e ~ja>Al e ~jco(i-l)At 


(19) 


The quantity e jmAl is constant for a given frequency a and constant sampling interval At . It follows that the 
discrete Fourier transform can be computed for a given frequency at each time step using one addition in Eq. (18) 

and two multiplications - one in Eq. (19) using the stored constant e J(0 1 for frequency a , and one in Eq. (18). 
There is no need to store the time-domain data in memory when computing the discrete Fourier transform in this 
way, because the data for each sample time is processed immediately. Time-domain data from the past can be used 
in all subsequent analysis by simply continuing the recursive calculation of the Fourier transform. In this sense, the 
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recursive Fourier transform acts as memory for the information in the data. More data from more maneuvers 
improves the quality of the data in the frequency domain without increasing memory requirements to store it. 
Furthermore, the Fourier transform is available at any time i At . The approximation to the finite Fourier transform 
is completed using Eq. (8). 

The recursive computation of the Fourier transform does not use a Fast Fourier Transform ( FFT ) algorithm 14 , 
and therefore would be comparatively slow, if the entire frequency band up to the Nyquist frequency \ 1(2 At) were 

of interest. However, rigid-body dynamics of aircraft lie in a rather narrow frequency band of approximately 
0.01-2.0 Hz. Since the frequency band is limited, it is efficient to compute the discrete Fourier transform using 
Eqs. (18) and (19) (which represents a recursive formulation of Eq. (7)) for only the selected frequencies 
a> k , k = 1,2 ,... ,M . With this approach, it is possible to select closely-spaced fixed frequencies for the Fourier 
transform and the subsequent modeling and still do the calculations efficiently. 

Using a limited frequency band for the Fourier transformation confines the data analysis to the frequency band 
where the dynamics lie, and automatically filters wide band measurement noise or structural responses outside the 
frequency band of interest. These automatic filtering features are important for real-time applications, where 
instrumentation error corrections and noise filtering would require additional computational resources. 

In past work on fighter aircraft short-period modeling, frequency spacing of 0.04 Hz on an interval of 
approximately [0.1-2] Hz was found to be adequate 9 n . Finer frequency spacing requires slightly more computation, 
but was found to have little effect on the results. When the frequency spacing is very coarse, there is a danger of 
omitting important frequency components, and this can lead to inaccurate parameter estimates. In general, a good 
rule of thumb is to use frequencies evenly spaced at 0.04 Hz over the bandwidth for the dynamic system. For good 
results, the bandwidth should be limited to the frequency range where the signal components in the frequency 
domain are at least twice the amplitude of the wide band noise level. However, the algorithm is robust to these 
design choices, so the selections to be made are not difficult. 

The recursive Fourier transform update need not be done for every sampled time point. Systematically skipping 
time points effectively lowers the sampling rate of the data prior to Fourier transformation. This saves computation, 
and does not have a significant adverse impact on the parameter estimation results until the Fourier transform update 
rate gets below approximately 5 times the highest frequency of interest for the dynamic system. The parameter 
estimation and covariance calculations in Eqs. (15)-(17) can be done at any time, but are usually done at 1 or 2 Hz, 
to save computations. Linearized aerodynamic characteristics rarely change faster than this, except in cases of 
strong nonlinearity, damage, failure, or rapid maneuvering. For these cases, the update rate can be increased, at the 
cost of more computations. 

Reference 6 explains that computing standard errors from the covariance matrix in Eq. (16) does not require 
correction for colored residuals. The standard errors computed from Eq. (16) are therefore a good representation of 
the error in the estimated parameters. Having high quality error measures is important for problems such as failure 
detection and control law reconfiguration. 

D. Data Forgetting 

The recursive Fourier transform in Eqs. (18) and (19) represents a data information memory for as long as the 
running sum is incremented. It follows that when the aircraft dynamics change, the older data should be discounted 
in some way, as has been done for time-domain approaches using a forgetting factor 6 . If this is not done, then the 
speed of response for the real-time parameter estimator is progressively degraded, as new information has to 
overwhelm an increasingly longer memory. Consequently, there is a trade-off between the desired rapid response of 
the parameter estimator to changes in the aircraft dynamics, versus retaining enough information from past data for 
sufficiently accurate model parameter estimates. 

If past values of the Fourier transform X t (®) computed from Eq. (18) are saved in computer memory, then it is 

possible to implement selective amnesia by simply subtracting past values of the running sum corresponding to the 
Fourier transform, or differences between past values of the running sum. For example, forgetting all data 
information content older than 10 sec (i.e., removing that data information content from the complex regression 
problem) could be implemented by subtracting the value of the running sums for the Fourier transforms at 10 sec 
ago from the current running sums. Similarly, to forget data information content collected between 5 and 7 sec ago, 
the difference between the running sums at 5 and 7 sec ago would be subtracted from the current running sum. The 
price to pay for this capability is the computer memory required to store past values of the running sums associated 
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with the Fourier transforms for each signal at each frequency. The memory requirements could be reduced by 
perhaps only saving the running sums at intervals of 0.5 sec, for example. 

The simplicity of Eq. (18) also makes it easy to see how exponential data forgetting could be implemented. In 
exponential data forgetting, each past value of the time-domain signal is multiplied by a forgetting factor A < 1 at 
each time step. In this way, old data is gradually devalued and eventually discarded. Usually A is chosen in the 
range 0.90 < A < 1.00 . To implement this, Eq. (18) is modified slightly to 


X t ( a ) = AX t _ x ( eo) + x(i ) e~ j(oiAt (20) 

and everything else remains the same as before. This simple approach is possible because the Fourier transform is 
linear with respect to the measured data x(i) . 

The challenge with data forgetting is not in the implementation, but rather in deciding how much data 
information content to forget, and when. At the present time, there are no concrete guidelines, so the choices are 
made based on analysis of results from simulation and flight data. 

E. Incorporating Prior Information 

Model parameter estimates can sometimes be improved by including information from modeling results based 
on prior data. In the case of real-time dynamic modeling, including prior information of this kind can reduce 
variations in the real-time parameter estimates, and improve convergence speed. 

One way to incorporate prior information is by using a mixed estimator formulation of the least squares cost 
function 6 . Assuming that the fit error variance cr p for the prior modeling is approximately equal to the fit error 
2 

variance <x for a model based on the current data alone, and denoting the vector of parameter estimates from a 
prior analysis by 9 p , with associated covariance matrix X p , the least squares cost function that incorporates this 

prior information is formulated as 6 


j(o)=-(z-xo)\z-xo) + -(o-o p ) T z- p x (o-o p ) 


The vector of parameter estimates that minimize this modified least squares cost function is 
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2 

where er is estimated from Eq. (17), as before. Reference 6 provides further details on the mixed estimator 
formulation of the least squares parameter estimation problem, and the associated solution. 

F. Optimized Excitation 

References 6, 11, and 16 describe a method for designing optimized inputs for use as control surface 
perturbations to excite aircraft dynamics. The form of each input is a sum of sinusoids with unique frequencies and 
phase shifts. Multiple inputs are designed to be mutually orthogonal in both the time domain and the frequency 
domain, and are optimized for maximum data information content in multiple axes over a short time period, while 
minimizing excursions from the nominal flight condition. The mutual orthogonality of the inputs allows 
simultaneous application of multiple inputs, which helps to minimize excitation time. Optimized inputs of this type 
were applied to the modified F-15B to collect data for real-time dynamic modeling. The optimized inputs were 
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applied as control surface excitations by summation with the actuator commands from the control system just before 
the actuator command rate and position limiting. Flight test examples are shown below in the Results section. 

III. Test Aircraft and Flight Data 


A. Fighter Aircraft Description 

The fighter aircraft used for this research is a pre-production Boeing F-15B that has been highly modified to 
support various test programs. The most visible modification is the inclusion of a set of canards near the pilot 
station, see Fig. 1 . The canards are a set of modified horizontal 
stabilators from a Boeing F/A-18 aircraft. The purpose of the 
canard addition was to increase maneuverability and load 
capability. An additional effect of the canards was to cause the 
aircraft to be statically unstable in the longitudinal axis at most 
subsonic speeds. The propulsion system consists of two Pratt 
& Whitney F100-PW-229 engines, each equipped with an 
axisymmetric thrust vectoring pitch/yaw balance beam nozzle. 

The thrust vectoring feature, however, was not used during the 
flight testing described here. Further information on the 
modified F-15B aircraft and associated flight test operations 
can be found in Ref. 17. 

Figure 1. Modified F-15B Fighter Jet Aircraft 

7. Control Surfaces 

The aircraft has five pairs of control surfaces: stabilators, canards, ailerons, trailing edge flaps, and rudders. 
Flaps and aileron droop are manually set by the pilot and only used for takeoff and landing configurations. 
Conventional pitch control is provided by symmetric deflection of the all-moving horizontal stabilators and canards. 
Roll control uses aileron and differential stabilator. Directional control is provided by rudder and differential canard 
deflection. Definitions of control surface deflections are given below. Trailing edge down is positive deflection for 
the wing and stabilator surfaces, and trailing edge left is positive for the rudder. 
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For the nominal flight control system, pilot stick and rudder inputs result in high correlation between the 
symmetric canard and angle of attack, rudder and differential canard, and differential stabilator and aileron. Flight 
data analysis showed that the symmetric canard and angle of attack exhibited nominal pair-wise correlation of 0.93, 
rudder and differential canard were correlated at 0.99, and the differential stabilator and aileron were perfectly 
correlated at 1.00. Given these high levels of correlated inputs, it would be impossible with pilot input maneuvers to 
distinguish, for example, between rolling moment generated by the aileron and rolling moment generated by the 
differential stabilator. Consequently, the complete stability and control derivative set cannot be obtained with pilot 
input maneuvers. 

2. Instrumentation and Data Acquisition 

The modified F-15B aircraft was equipped with a research-quality instrumentation system. A nose boom was 
installed and calibrated for free stream pitot-static and flow angle measurements. An inertial instrumentation 
package provided 3-axis linear accelerometer and angular rate measurements. Heading, pitch, and bank angle 
measurements were obtained from the aircraft inertial navigation system. Control surface positions were measured 
using variable differential transformer sensors. Fuel measurements were obtained from the three fuselage and two 
wing tanks. Fuel measurements were used to compute total aircraft weight, center of gravity, and mass moments of 
inertia. Mass and geometry characteristics of the aircraft are given in Table 1. Sensor positions and the center of 
gravity (c.g.) location were used to correct the linear accelerometer and flow angle measurements to the c.g. 
Angular accelerations were not measured. Flight data was collected at 40 Hz. 
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B. S-2 Subscale Jet Transport Aircraft Description 

The S-2 aircraft is a subscale model of a jet transport 
aircraft A photograph of the aircraft in flight is shown in 
Fig. 2. The subscale aircraft has a single jet engine mounted in 
the aft fuselage and retractable tricycle landing gear. Mass and 
geometry characteristics of the aircraft are given in Table 1. 
Further information on the S-2 subscale jet aircraft and 
associated flight test operations can be found in Ref. 18. 

1 . Control Surfaces 

Control surfaces on the aircraft are conventional elevator, 
aileron, rudder, and inboard trailing-edge flaps. Definitions of 
control surface deflections are given below. Trailing edge 
down is positive deflection for the wing and elevator surfaces, 
and trailing edge left is positive for the rudder: 


$ 



Figure 2. S-2 Subscale Jet Transport Aircraft 


S=-(s„ +S e ) S=-(s a -8 a ) 
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(25) 


The aircraft can be flown by a safety pilot using direct visual contact and conventional radio control. A research 
pilot executes flight test maneuvers from inside a mobile control room, using a synthetic vision display drawn from 
telemetry data and a local terrain database. Inputs from the research pilot and ground-based flight control are 
transmitted to the aircraft by telemetry to the aircraft control surface servomechanisms. 

2. Instrumentation and Data Acquisition 

The S-2 aircraft was equipped with a micro-INS, which provided 3-axis linear accelerometer measurements, 
angular rate measurements, estimated attitude angles, and GPS velocity and position. Air data probes on each 
wingtip (visible in Fig. 2) measured angle of attack, sideslip angle, static pressure, and dynamic pressure. 
Measurements from static pressure sensors and ambient temperature sensors were used to compute air density and 
altitude. Engine speed in rpm was measured and used as input to an engine model to compute thrust. The engine 
model was identified from ground test data, with adjustments for ram drag identified from flight data. 
Potentiometers on the rotation axes of all control surfaces measured control surface deflections. Mass properties 
were computed based on measured fuel flow, pre-flight weight and balance, and careful inertia measurements of the 
aircraft on the ground. The pilot stick and rudder commands and throttle position were also measured and recorded. 
Flight data was collected at 50 Hz. 


IV. Results 

A. Data Information Requirements for Real-Time Dynamic Modeling 

In this section, results will be presented from investigations regarding data information content necessary for 
accurate real-time modeling. 

Figure 3 shows a relationship between signal-to-noise ratio of aircraft measured responses to the mean error in 
real-time parameter estimates computed using the FTR method described earlier. The data for this analysis was 
generated with a linear simulation of the modified F-15B lateral dynamics at a flight condition of Mach 0.75, trim 
angle of attack 2 deg, and 20,000 ft altitude. Optimized multi-sine inputs were applied simultaneously to the lateral 
control surfaces, as shown in Fig. 4. Repeated simulated data runs with different output signal-to-noise ratios were 
generated by uniformly reducing the amplitudes of the inputs shown in Fig. 4, generating new simulated outputs, 
then adding a single realization of noise sequences extracted from flight data to the simulated outputs. This caused 
the output signal-to-noise ratios to vary in a uniformly decreasing fashion. For each data run, the FTR method was 
applied, with no prior information, to produce real-time parameter estimates. Final parameter estimates and 
standard errors from each 18 sec run were used to generate the data points for Fig. 3. Because the problem involves 
multiple parameters and multiple measured outputs, average values of the parameter estimate accuracies and the 
output signal-to-noise ratios were used to make the simple summary plot shown in Fig. 3. Based on Fig. 3, 
signal-to-noise ratios of the measured aircraft responses should preferably be 10 or more, and not less than about 3, 
for reasonable accuracy of the real-time parameter estimates. This finding is important in designing excitation 
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inputs, and also suggests that lower levels of excitation could be used to obtain accurate real-time dynamic modeling 
if noise levels were decreased. It is of interest to obtain good real-time parameter estimation results with the lowest 
possible signal-to-noise ratio on the measured outputs, so that the aircraft and pilot are disturbed as little as possible 

by any excitation applied to the control surfaces. For 
output signal-to-noise ratio equal to 10, the input 
amplitudes were 0.24 deg, and the mean parameter 
estimate error was 2.7 percent. This choice of input 
amplitudes produced excellent parameter estimate 
accuracy, while keeping the root-mean-square value of 
the lateral acceleration a y below 0.02 g. This 

magnitude for the root-mean-square lateral 
acceleration is less than the typical value for moderate 
turbulence, which is 0.05 g. 

Similar optimized multi-sine inputs shown in 
Fig. 5 were applied to the modified F-15B aircraft in 
flight at the same flight condition, Mach 0.75, 2 deg 

trim angle of attack, and 20,000 ft altitude. The 
Output Signal-to-Noise Ratio aircraft is i ongitudina u y unstable at this flight 

Figure 3. Dependence of mean real-time parameter condition, so feedback control was operating when the 
estimation error on output signal-to-noise ratio control surface excitations were applied by summation 

with the actuator commands from the control system. 
As shown in Fig. 5, all longitudinal and lateral control surfaces were moved simultaneously. The dashed lines show 
the optimized multi- sine input design, where the inputs are mutually orthogonal in both the time and frequency 
domains. The solid lines show the actual control surface measurements from flight, which are distorted from the 
optimal waveforms. The highest pair-wise correlation for control surface deflections is between stabilator S s and 
canard S c , because of the high-gain feedback necessary to compensate for the longitudinal instability at this flight 

condition. Nevertheless, the control surface deflections were largely uncorrelated, as shown by the pair-wise 
correlation coefficients listed in Table 2. These correlation levels are well below the 0.9 correlation level typically 
used as an upper limit for allowable correlation of explanatory variables 6 . Although the relatively high-gain 
feedback control ruined the input orthogonality, the data still exhibited very low correlations among the control 
surface deflections, which is desirable for accurate dynamic modeling 6 . 

Using optimized multi-sine inputs, the effect of input amplitude was investigated by repeating maneuvers with 
varying input amplitudes on the same flight. Figure 6 shows the variation in amplitudes for the lateral control 
surface deflections in flight. Selected real-time parameter estimates are shown in Fig. 7. The markers (X,0,+) 
indicate the real-time parameter estimates, while the lines only connect these markers to more clearly show the 
evolution of the parameter estimates. Real-time modeling was done using the FTR method with no prior 
information, frequency vector / = [0.10 0.12 ... 2.00] Flz , recursive Fourier transform run at the data rate of 

40 Flz, and real-time parameter estimation calculations were done at 2 Hz. Selections for the frequency vector and 
update rates were not critical to the results produced - similar results could be obtained for a fairly wide-range of 
choices of these quantities. The results shown in Fig. 7 are typical of all of the real-time parameter estimation 
results. The excitation amplitudes made little difference in either the speed of convergence or the final parameter 
estimates. In general, control derivative estimates (e.g., C t ) were the fastest to converge and usually the most 

accurate, followed by static stability derivative estimates (e.g., C„^ ) and damping parameter estimates (e.g., Cj ). 

Roughly 8 sec of excitation was required for all model parameters to converge, with no prior information. Standard 
errors are not shown in Fig. 7, to reduce clutter on the plots. 

The mean signal-to-noise ratios for the measured lateral aircraft responses to low, medium, and high amplitude 
inputs were 6.6, 8.5, and 11.4, respectively, with mean parameter estimate errors 15.9, 12.0, and 6.4 percent, 
respectively. These results have the same character as the data shown in Fig. 3, but with mean parameter errors 
roughly 3 to 4 times higher. This is not unexpected, because flight data has noise components from structural 
responses, nonlinearities, and model structure errors that are not present in the linear simulation. This leads to 
higher mean parameter errors for flight data analysis and modeling. The root-mean-square values of the lateral 
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acceleration a y for low, medium, and high amplitude inputs were 0.020 g, 0.028 g, and 0.047 g, respectively, 
compared to 0.05 g, which is typical for moderate turbulence. 

Variation in the frequency content of the excitation inputs was studied by applying two different optimized 
excitation input designs in flight, one for a short (6 sec) time length and consequently more sparse frequency 
content, and one with a longer (15 sec) time length and richer frequency content. A comparison of the frequency 
content for the two optimized input designs is shown in Fig. 8. Both optimized input designs were flown on the 
modified F-15B aircraft, at the same flight condition as before. Typical real-time parameter estimation results are 
shown in Fig. 9, with 95 percent confidence intervals based on estimated parameter standard errors shown by the 
vertical bars through each estimate. In general, the more sparse frequency content led to slightly poorer 

convergence and accuracy, except for control derivatives like C n . Figure 9 shows that a more effective approach 

for short excitation times is to truncate a maneuver designed for a longer time length, because the longer maneuver 
has a wider variety of frequencies, more closely spaced along the frequency axis, cf. Fig. 8. 

Figure 10 shows a comparison of selected real-time parameter estimates based on flight data from the modified 
F-15B aircraft at the same flight condition. Optimized multi-sine excitation was done first by moving only the 
longitudinal control surfaces (<5 V and 8 C ) simultaneously, then again moving all longitudinal and lateral control 

surfaces (8 S , S c ,S a ,S r ,8 dc , and simultaneously. The results in Fig. 10 show that there was essentially no 

difference in the real-time parameter estimation behavior for these cases, indicating that the multiple input design 
was effective, the equation-error parameter estimation worked as expected on each equation individually, and finally 
that there was negligible aerodynamic coupling for the maneuvers resulting from these optimized excitation inputs. 
Similar results were found when comparing lateral excitation alone to combined longitudinal and lateral excitation. 
Real-time parameter estimate behavior for other model parameters was similar to what is shown in Fig. 10 for the 
three selected longitudinal parameters. 

The results in Fig. 10 indicate that combined longitudinal and lateral maneuvers with all control surfaces moving 
at once should be used for maximum efficiency in real-time dynamic modeling. This approach can provide a full set 
of stability and control modeling results in a single relatively short maneuver. The mutual orthogonality of the 
multiple input design makes this possible. Finally, Fig. 10 shows that the FTR method has no trouble identifying 

dynamic model parameters when the aircraft is statically unstable longitudinally (C m > O) , as in this case. 

B. Data Information Management 

The real-time dynamic modeling results shown in Figs. 7, 9, and 10 were generated with the FTR algorithm 
assuming no priori information. If instead the FTR algorithm is provided some information on the values of the 
model parameters, in terms of parameter estimates and associated covariance matrix from a prior analysis, the 
convergence behavior of the real-time parameter estimates from the FTR method changes. Figure 1 1 shows real- 
time parameter estimate histories for a combined longitudinal and lateral maneuver using medium-amplitude 
optimized multi-sine excitation and no prior information, along with real-time parameter estimate histories using the 
same data, but with prior information from parameter estimates and associated covariance matrix at the end of a 
low-amplitude combined maneuver. Real-time estimate variations are greatly reduced with the addition of prior 
information, particularly at the beginning of the maneuver. However, this regularization also results in something 
like a weighted average for the final parameter estimate, which is different from what would be computed with no 
prior information, as can be seen clearly in the top plot of Fig. 1 1 for C m . The extent of this effect depends on the 

values and accuracy of the parameter estimates provided as prior information. 

The quality of the estimated parameter error bounds was judged by how well the error bounds characterized the 
difference between the associated real-time parameter estimate and the final converged parameter estimate. The 
computed error bounds accurately characterized these discrepancies, which is consistent with similar findings in 
previous work 6 ' 13 . Although such a measure is easy to compute post-flight, metrics that can be computed in real 
time must be devised. Adequacy of local real-time dynamic models can be assessed using short-term prediction 
tests or tests involving the size of the estimated parameter standard errors. Although various metrics have been 
proposed, research continues toward devising a generally-applicable metric for real-time model validation. 

C. Real-Time Dynamic Modeling in Operational Flight 

All of the flight testing on the modified F-15B aircraft was at a single flight condition, using various optimized 
excitation inputs. To investigate the performance of the real-time dynamic modeling method in operational flight 
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with low data information content, flight data from the S-2 subscale jet transport aircraft in normal approach (i.e., 
without specific excitation inputs) was used. Figure 12 shows real-time parameter estimation results for the S-2 
aircraft during approach in landing configuration (full flaps, landing gear down). The real-time parameter estimates 
shown in Fig. 12 were obtained with no prior information. 

The real-time dynamic modeling results in this case are of similar quality to the results obtained with optimized 
multi-sine excitation inputs. Part of the reason for this is that the S-2 aircraft must be actively controlled by the pilot 
during approach to landing, and there is no feedback control system operating. However, these results show 
encouraging success in estimating real-time dynamic model parameters without specific control surface excitation. 
Further flight research is planned for the S-2 aircraft, which has been equipped with the capability to apply arbitrary 
excitations to any or all control surfaces. This capability will be used to investigate lower limits on signal-to-noise 
ratios for real-time dynamic modeling in changing flight conditions and for real-time dynamic modeling near stall 
and departure. 

Real-time dynamic modeling was also applied in S-2 flight testing during 3-axis pilot doublet maneuvers (pitch, 
yaw, and roll doublets, in sequence), with results displayed inside the mobile operations station (MOS) 18 in real 
time. Figure 13 shows the real-time display for selected parameter estimates during a 3-axis pilot doublet maneuver 
at 80 knots and 800 ft altitude. The dashed lines (added post-flight) indicate values of the parameter estimates 
obtained using a time-domain output-error method that employs an iterative nonlinear optimizer in post-flight batch 
mode 6 . It is evident that after only a few seconds of data, the real-time dynamic modeling algorithm (operating with 
no prior information) computed very accurate stability and control derivative estimates. Note that each stability or 
control derivative is inaccurate until the relevant measured quantity shows some significant variation (i.e., data 
information content), after which the associated derivative estimate converges quickly. Finally, the 95 percent 
confidence intervals for the parameter estimates, indicated by the vertical bars through each parameter estimate 
marker, properly characterize the accuracy of the real-time parameter estimates. 

V. Concluding Remarks 

Simulated data and flight test results were used to investigate data information requirements for accurate real- 
time dynamic modeling. This technology provides the capability to identify local dynamic models in real time from 
noisy data records. With proper data information management, the approach could be applied to changing flight 
conditions that occur throughout a flight, and/or could identify local dynamics for cases of failure or damage that 
can be characterized by model parameter changes. 

A real-time dynamic modeling algorithm was described in detail. The algorithm uses equation-error in the 
frequency domain, and a recursive Fourier transform. Modifications necessary to incorporate information from a 
previous analysis into the least squares formulation, or to implement data forgetting in the frequency domain, were 
developed and explained. 

Results from flight investigations on a modified F-15B aircraft showed that optimized multi-sine inputs provided 
excellent data information content for real-time dynamic modeling, even with a relatively high-gain feedback 
control system operating. Mean signal-to-noise ratio of aircraft flight responses was inversely correlated with mean 
real-time parameter estimation error, based on analysis of both simulation data and flight data from the modified 
F-15B aircraft. This suggests that smaller excitations can be used as measurement noise levels get smaller, making 
the required excitations less noticeable to the pilot. Refinements to the optimized excitation input designs were 
identified based on flight results. Combined longitudinal and lateral optimal excitations were found to efficiently 
provide data information content for real-time dynamic modeling in all axes simultaneously. Incorporating 
information about the dynamic model parameters from a prior analysis into the real-time modeling algorithm was 
shown to effectively reduce variability in real-time parameter estimates, particularly at the start of a maneuver. 
Real-time estimates of parameter errors computed from flight data were found to properly characterize the accuracy 
of the real-time parameter estimates. Finally, flight data from a subscale jet transport aircraft was used to show that 
the real-time dynamic modeling can also work well for 3-axis pilot doublet maneuvers, as well as ordinary flight 
operations maneuvers, such as approach to landing. 

Accurate local real-time dynamic modeling results can be used as building blocks for global real-time modeling. 
Other uses for the local real-time dynamic modeling results are in improving flight test efficiency, and prediction of 
dynamic behavior and departure. Identified local linear models can be also used in real time to track eigenvalues 
and eigenvectors of the dynamic system. This information can be used in bifurcation analysis for nonlinear flight 
regimes. 
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Table 1 Geometry and mass properties 



Modified 

S-2 


F-15B 

length c , ft 

15.94 

0.908 

wing span b , ft 

42.70 

7.083 

wing area S , ft 2 

608.0 

7.046 

x ref , in 

557.2 

42.628 

Jref > in 

0.000 

0.000 

z re f > in 

116.3 

0.000 

x cg , in 

560.40 

42.728 

ycg ’ in 

0.14 

0.000 

z c> , , in 

117.41 

0.519 

m , slugs 

1234 

1.502 

I x , slugs-ft 2 

24,830 

1.077 

I , slugs-ft 2 

196,225 

4.163 

I, , slugs-ft 2 

216,155 

5.016 

I xz , slugs-ft 2 

-5329 

0.416 


Table 2 Pair-wise correlation matrix for control surface deflections 



8 S 


8 a 

S r 

Sdc 

Sds 


1 

-0.293 

-0.085 

-0.026 

0.126 

0.016 

8 C 


1 

-0.120 

-0.065 

-0.043 

-0.149 

s a 



1 

0.083 

0.011 

-0.010 

5 r 




1 

0.177 

-0.029 

^dc 





1 

-0.059 


$ds 1 
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Figure 4. Simulated data for optimized multi-sine inputs, modified F-15B aircraft 
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Figure 5. Optimized multi-sine input distortion, modified F-15B aircraft 
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Figure 6. Input amplitude variation, modified F-15B aircraft 



Figure 7 Real-time parameter estimates for varying input amplitudes, 
modified F-15B aircraft 
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Figure 8. Excitation input spectra for 6 sec and 15 sec longitudinal maneuvers, modified F-15B aircraft 





Figure 9. Real-time parameter estimates for varying frequency content, 
modified F-15B aircraft 
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Figure 10. Real-time parameter estimates for combined 
maneuver compared to longitudinal only, modified F-15B aircraft 
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Figure 11. Real-time parameter estimates with and without 
prior information, modified F-15B aircraft 
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Figure 12. Real-time parameter estimates during approach in 
landing configuration, S-2 aircraft 
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Figure 13. Real-time parameter estimates for a 3-axis pilot doublet maneuver, S-2 aircraft 
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